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Abstract 

We describe a non-perturbative procedure for solving from first principles the light-front 
Hamiltonian problem of SU{N) pure gauge theory in D spacetime dimensions (-D > 3), based on 
enforcing Lorentz covariance of observables. A transverse lattice regulator and colour-dielectric 
link fields are employed, together with an associated effective potential. We argue that the light- 
front vacuum is necessarily trivial for large enough lattice spacing, and clarify why this leads to 
an Eguchi-Kawai dimensional reduction of observables to 1 + 1-dimensions in the N —> oo limit. 
The procedure is then tested by explicit calculations for 2 + 1-dimensional SU (oo) gauge theory, 
within a first approximation to the lattice effective potential. We identify a scaling trajectory 
which produces Lorentz covariant behaviour for the lightest glueballs. The predicted masses, 
in units of the measured string tension, are in agreement with recent results from conventional 
Euclidean lattice simulations. In addition, we obtain the potential between heavy sources and 
the structure of the glueballs from their light-front wavefunctions. Finally, we briefly discuss 
the extension of these calculations to 3 -I- 1 dimensions. 



^On leave from DAMTP, University of Cambridge, Cambridge CBS ODS, England. 
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1 Introduction 



Light-front Hamiltonian quantisation is a very promising approach to the non-perturbative calcu- 
lation and understanding of hadronic substructure, as a function of the quark and gluonic degrees 
of freedom. This is both of intrinsic interest and important for the proper identification of new 
physics beyond the Standard Model. Recently, we have developed an idea originally suggested 
by Bardeen and Pearson many years ago ^, for formulating the light-front Hamiltonian of QCD 
on a transverse lattice. The continuum limit is taken in time and one space direction, while the re- 
maining transverse degrees of freedom are represented as colour-dielectric link variables in a gauge 
invariant way on the transverse lattice. Since light-front wavefunctions carry the non-perturbative, 
coherent information about hadronic bound and scattering states in a general Lorentz frame, a so- 
lution to this problem would open up a vast amounts of existing and future experimental results to 
theoretical scrutiny ||]. We believe that the menagerie of light-front Hamiltonian, transverse lattice, 
and effective colour-dielectric fields, is a very potent combination in the analysis of the strong inter- 
actions at intermediate scales. Each contributes to the solution of long-standing technical problems 
associated with the others. In addition, profound simplifications take place when one considers only 
the leading order of the expansion in the number of colours N. In this paper, we describe 
in detail how non-perturbative first principles calculations may be performed for pure gauge the- 
ory, in 2 -|- 1 and 3-1-1 dimensions, by imposing Lorentz covariance on physical observables. The 
procedure is then tested with explicit numerical and analytic computations for 2 -f- 1 dimensions at 
large iV, comparing successfully with recent 'state-of-the-art' results on the glueball spectrum from 
conventional Euclidean lattice simulations, as well as producing many new results. 

The choice of dynamical variables and regulator in canonical light-front quantisation is crucial 
to a practical solution. The well-known Lorentz boost invariance of light-front wavefunctions comes 
with a price, because, for example, the light-front quantisation surface is not rotationally invariant; 
typically the regulators chosen will break rotational invariance. An effective potential must be 
introduced to restore full Lorentz invariance. This effective potential plays another important role 
to encode what is usually, id est in equal-time quantisation, termed 'vacuum structure.' One of the 
attractive features of light-front quantisation is that, in the true groundstate of a theory, the light- 
front vacuum is trivial — there can be no pair creation — except possibly for some massless gauge 
modes of zero momentum. But in the true groundstate, much of the complication usually attributed 
to the vacuum has only been shifted into effective interactions appearing in the Hamiltonian. 

Typically a light-front Hamiltonian will consist of a 'bare' part derived from the classical 
action, and the effective potential consisting of countcrtcrms He- The role of the latter is to enforce 
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the fundamental (unbroken) symmetries required of the theory, which are broken by the choice of 
regulators. For practical purposes, one must decide on some physical grounds which subset of all 
possible operators one is going to consider in He , a subset which will then be systematically enlarged 
as greater accuracy is demanded. The method we will use in our work, to fix the couplings to this 
subset of operators, is to solve the entire theory as a function of the couplings, which are then tuned 
by demanding unbroken symmetries of calculated observables. This straightforward approach will 
allow direct calculations involving typical hadronic scales, intermediate between the asymptotically 
free region and nuclear scales. Naively, one might think that the approach is computationally 
unviable, but this is not the case if one uses Hamiltonian methods. There already exists at least 
one successful application of this approach to fixing counterterms. Parity, which is not a manifest 
symmetry of the light-front, is typically broken by regulators. Burkardt Q] has shown how, in the 
case of the Yukawa model in 1 + 1 dimensions, a simple counterterm coupling can be adjusted to 
restore parity invariance in a suitable observable. The general problem in higher dimensions involves 
the far more complicated task of ensuring the continuous Lorentz symmetries also. This will form 
the main part of our job in solving pure gauge theory. 

To realise the above scenario, we apply light-front quantisation to a set of massive effective degrees 
of freedom on the transverse lattice — colour-dielectric link variables — which efficiently organise 
a gauge invariant Hamiltonian in the true vacuum. These variables are supposed to represent an 
average over the high frequency transverse gluon degrees of freedom. The fact that (lattice) gauge 
invariance is not broken helps to restrict the number of counterterms needed in He, and leads to 
a simple understanding of confinement. The dielectric variables will provide the required physical 
basis for systematic approximation to the low-energy effective potential. Moreover, they are flexible 
enough that in the light-front Hamiltonian formulation one can solve for physical observables, with 
reasonable computational effort, as a function of the couplings in the effective potential. This 
formulation therefore satisfies the requirements necessary for our approach. In the first half of this 
paper, we construct in detail the framework for making calculations in pure gauge theory at large 
N. This framework is developed in 3 + 1 dimensions, but can also be applied in 2 + 1 dimensions 
with some trivial modifications. 

As a quantitative test in a highly relativistic, strongly interacting system, we then perform 
detailed calculations in 2 + 1 dimensions. This choice as a testing ground is motivated by a number 
of reasons. A more reliable computation can be performed in 2 + 1 dimensions than 3 + f dimensions. 
Teper has recently completed extensive and precise conventional Euclidean lattice Monte Carlo 
(ELMC) simulations in three dimensions for various N. The resulting spectra may be extrapolated 
to large N, providing a benchmark for our work. Finally, the physics of Yang-Mills theory in 2 + 1 
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dimensions is remarkably similar that in 3 + 1 dimensions. In our previous work ^ (hereafter referred 
to as I), we studied the effective potential in 2 + 1 dimensions by matching our glueball spectrum 
to Teper's preliminary data. For a truncated effective potential, we were able to identify a coupling 
constant trajectory along which mass ratios were approximately constant and agreed reasonably well 
with the ELMC continuum predictions. This is a perfectly acceptable phenomenological approach, 
but one would obviously like to avoid fitting to an existing spectrum, whether taken from experiment 
or another theoretical method. The first principles calculations we perform in this paper use only 
general concepts like gauge and Lorentz invariance to fix the effective potential. In this way, we 
can make independent predictions of all physical quantities. The contents of the present paper are 
summarised below. 

In the next section, we briefly review the colour-dielectric formulation of lattice gauge theory and 
its light-front Hamiltonian limit in order to set notation. In Section ^ we carry out the light-front 
quantisation, paying particular attention to the interplay of cutoffs, gauge fixing, and energetics. 
This allows us to show that in the 'colour-dielectric regime' where the light-front Hamiltonian is 
quantised, light-front zero modes may be neglected, id est the light-front vacuum is essentially trivial 
(this justifies the naive light-front gauge choice often used). We observe that linear confinement is 
generic to the theory in the colour-dielectric regime. All our work is performed in the large- limit 
because of remarkable simplifications which occur. In Rcf. § and I it was proved that, in this 
limit, one effectively solves the light- front Hamiltonian of a I + I-dimensional gauge theory coupled 
to 2^^^ adjoint matter fields, for D space-time dimensions. We clarify in Section |^ the physical 
role of this Eguchi-Kawai dimensional reduction. Unlike lattice path integral quantisation, where in 
the leading order in N vacuum structure undergoes dimensional reduction to zero dimensions, the 
leading order in N of the transverse lattice light-front quantisation exhibits dimensional reduction 
of physical observables to I + I dimensions. Section ^ outlines procedures for measuring the string 
tension in the transverse lattice direction via winding modes, and the full spatial potential between 
two heavy sources, extending work of Ref. [lO[ |. In Section we construct states with nonzero 
transverse momentum. We use the resulting dispersion relation, together with the string tension 
measurements, to construct a simple algorithm for fixing the effective potential based on Lorentz 
covariance. 

We then turn to a numerical implementation of our procedure for 2 + 1-dimensional gauge theory 
(Section 0). Most of the relevant formulae and descriptions of our numerical methods have been given 



in I; we review and extend them in Sections B.l, 8.2, and the Appendices. Section ^ presents our 
results for the scaling behaviour, glueball masses and wavefunctions, the string tensions and heavy- 
source potential. Finally, we conclude with a brief discussion of the extension to 3 + 1 dimensions 
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of these calculations. 

2 Colour-Dielectric Lattice Gauge Theory. 
2.1 Euclidean Lattice 

For the purposes of studying a light-front Hamiltonian limit, Wilson's formulation of lattice gauge 
theory leads to difficulties due to the non-linearity of the link variables We therefore take 
as our starting point Mack's linearised 'colour-dielectric' formulation of lattice gauge theory Q. A 
similar lattice gauge theory was first introduced by Weingarten as a generator of the spacetime 
lattice regularisation of the Nambu string action. 

On a 4-dimensional hypercubic lattice of spacing a, introduce a complex matrix M^{x) € 
GL{N,C) for the link attached to site x = {x'^,x^,x^,x^} in the direction x^^. The partition 
function is 




with colour indices i,j G {1, • • • , N} and M^{x) — M_^{x + ajl). Throughout, fl is the unit vector 
in the x^ direction and — /i = — /i. For an appropriate choice of action 5, this theory is invariant 
under lattice gauge transformations 

M^,{x) U\x)M^,{x)U{x + ajl) , U £ SU{N) . (2) 

Since we will eventually be interested in an expansion of the action in powers of the field M in the 
large N limit, we take S[M] to consist of Wilson loops and their products. Note that this will include 
Wilson lines which 'backtrack', such as Tr{ Af^t(2;)Mt (a;)} , because M is not necessarily a unitary 
matrix. (Determinants of M are also gauge invariant, but involve an infinite product of fields when 
N = oo.) Since we are interested in QCD, S is tuned in such a way that, as the correlation length 
approaches infinity in direction x^ (lattice spacing a — + in physical units) , 

\[m MJx) = 1 + iaAf, + Oia"^) , (3) 

a—*0 

where A^^ is the usual hermitian gauge potential. Thus, M must tend asymptotically to the unitary 
Wilson link variable and in the isotropic continuum limit S will be proportional to Trji^^t^F^''}. 
For future reference, it is convenient to introduce lattice versions of the covariant derivative and field 
strength tensor 

Df,M^{x) = M^{x)NU{x + op) - AU{x)M^{x + dV) . (4) 
This natural definition reduces to -F^^ in the continuum. 
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Figure 1: All Wilson loops and their local products of total length < 4. To those shown above one 
must add the different orientations of each loop in D spacetime dimensions and reversal of arrows 
(only the plaquette gives a distinct loop in the latter case.) In the case of products of Wilson loops, 
we have only shown the case where the loops lie exactly on top of one another (local product). 



More generally the coefficients of S are to be chosen to obtain scaling, id est Green's function 
observables are to be independent of the correlation length(s). We will demonstrate existence^ 
of an approximate scaling trajectory, in the light-front Hamiltonian limit, by demanding Lorentz 
covariance of boundstate solutions to the theory with S truncated to a finite number of terms. This 



empirical approach is similar in spirit to the Monte Carlo Renormalisation Group 1 14 . However, 
we find it more convenient to directly minimise violations of Lorentz covariance in observables, as a 
function of couplings in S, rather than try to perform explicit 'block-spin' transformations. 

The Wilson loops up to length 4 are illustrated in Fig. |^. The quadratic loop Tr\^Mf^{x)Mj^{x)} , 
id est the 'mass squared' term of the M^{x) fields, evidently plays a special role. It must obviously 
have negative (tachyonic) coefficient in the continuum limit (^). If it has a positive coefficient, a 
'strong coupling' expansion of the path integral ([T]) about M^{x) = may be performed |jl^. If 

^Uniqueness would follow from the assumption that the continuum limit of QCD is unique, since the light-front 
Hamiltonian approach requires us to take the continuum limit (^) in two directions and then demand approximate 
Lorentz invariance. 



the scaling trajectory passes through such a region, Mack showed that a colour-dielectric picture of 
confinement results Accordingly, we shall refer to this as the colour- dielectric regime. There is 

evidence for SU{2) that a partition function of the form (|l]), subject to 'block-spin' transformations, 
does have renormalisation group trajectories which pass into the colour-dielectric regime at short 
enough correlation length p6[ |. Therefore, the picture one should keep in mind is the following. If one 
decomposes M^{x) — HU into a hermitian matrix H and a unitary matrix U , in the continuum limit 
a ^ 0, [/ is identified with Wilson's link variable. H = Hq + H gets a VEV Hq, while the fluctuation 
H becomes very heavy and decouples. Near the continuum limit, Hq appears in the equations of 
motion like a generalised dielectric constant jj]. In this regime the field M^{x) is tachyonic. As the 
lattice spacing a is increased, a scaling trajectory in S pushes one into a region of positive mass 
squared for M^,(x), where iJo vanishes and H is fully dynamical. The mass of M^{x) then increases 
with a. We will find further evidence for this in our light-front Hamiltonian calculations later, where 
it is rather crucial to work in a regime of positive mass squared for M^{x). 

2.2 Light-Front Continuum Limit. 

It is difficult to perform quantitative calculations with the the dielectric lattice gauge theory de- 
scribed above because one must search the multi-dimensional coupling constant space of S for a 
scaling trajectory. We will show that these problems can be overcome in the light-front Hamiltonian 
limit. To employ a light-front Hamiltonian formulation, it is desirable to take the continuum limit in 
two directions. Thus, in D spacetime dimensions we will take the lattice spacing to zero in the x'^ and 
x^~^ directions while keeping the lattice spacing a in the D — 2 remaining "transverse" directions 
constant. Throughout, boldface will denote vectors in the transverse directions, b, with components 
where r,s E {1, . . . , D — 2}. Continuum Lorentz indices are denoted by a, /3 G {0, D — 1}; thus 
we write the inner product of two Lorentz vectors as a'^bp, = a°'ba — a • b. Transverse coordinates x 
and y will hereafter be restricted to values corresponding to transverse lattice sites. The covariant 
derivative in direction x°' becomes 

UaMr(x) = {da + «^a(x)) A/^ (x) - zAfr(x)Aa(x + flf) , (5) 

in the partial continuum limit. One finds a transverse lattice action 

j dx'^dx^ (Tr|;D„Af,(x) (;d"M,(x))^| - ^ Tr{FapF^f'} ~ V^[M]^ , (6) 

where we have chosen to normalise A/r(x) to have a canonical kinetic term, and repeated indices 
are summed unless states otherwise. The dimensionful coupling G^{a) is such that a^^^G^ gjj 
in the classical continuum limit a — > 0, where gj^ is the continuum gauge coupling in D dimensions. 
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If we keep only Wilson loops and their products of total length 2 and 4 (which we justify in the next 
section) , the effective potential generically consists of Wilson loops among those shown in Fig. |l|, 
but now understood to be in the remaining D — 2 dimensional transverse lattice and pinned at one 
point to X. 

The basis for an effective QCD light- front Hamiltonian was earlier proposed in Ref. Ql, 
motivated by considerations of two-dimensional sigma models. Subsequently, an exploratory glueball 
calculation was performed Q for some simple choices of V . Those authors speculated that one might 
be able to find a scaling trajectory in V, corresponding to QCD. We have extended the methods 
of Ref. 1^ in order to answer this question, and perform in this paper quantitative first principles 
calculations. This allows us to confirm the existence of an approximate scaling trajectory at large 
transverse lattice spacing (at least in 2 -I- 1 dimensions), and show that it yields the same results for 
physical observables that are available from conventional methods of studying QCD. 

As mentioned above, to take best advantage of the simplifications on the light-front afforded by 
linear link variables Mr(x), it is quite important to be able to quantise about Mrijn) = 0. Our 
calculations therefore must implicitly test a light-front Hamiltonian version of the colour-dielectric 
regime: 

In the light-front Hamiltonian limit of lattice action S , we search for a scaling trajectory by canonical 
quantisation about Mr{x) = 0, taken as the true groundstate. 

Once outside of the purported colour-dielectric regime, when a is sufficiently small and Mr(x) = 
is no longer the minimum of energy, we expect the zero mode structure of light-front quantisation 



to become much more complicated [17 



3 Canonical Light-Front Quantisation. 

3.1 High Energy Cutoffs 

We introduce light-front co-ordinates = ~ {x^ ± a;'^^^)/\/2, for all two- vectors in the 
{x^ , x^~^) plane, and quantise by treating x'^ as canonical time. The variables conjugate to x^ 
and X are = (fc° =F /c^~^)/\/2 and k respectively. Particles and antiparticles have /c+ > 0. The 
transverse lattice theory will have divergences associated with the /c+ — > limit. This is a high k^ 
light-front energy limit, as follows from the dispersion relation 

for free particles of mass fi. At finite transverse lattice spacing a, this is the only source of divergence, 
since k is bounded by the Brillouin zone ir/a. In fact, the theory is just a set of coupled 1-1-1 
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dimensional continuum field theories, and there are only the mild free-field divergences characteristic 
of two-dimensional super-renormalisable theories. For pure gauge theory, finiteness is ensured by a 
suitable normal-ordering prescription. 

At the level of free fields, it follows that if /i^ > the quantum vacuum is free of particles. Because 
light-front momentum A:+ > and is conserved, the zero momentum vacuum state with even one 
massive (fc+ = 0) particle present is infinitely higher in energy than one with no particles. In the 
interacting theory, the particles must not condense either. In other words, one assumes that one is 
already expanding about the true vacuum state, and all condensation effects are parameterised by 
effective interactions in the Hamiltonian. This scenario is precisely realised in the colour-dielectric 
regime; the dielectric fields Afr(x) are massive and A/r(x) = is the true vacuum. The non- 
trivial problem is to demonstrate a coupling trajectory in this regime which restores all fundamental 
(unbroken) symmetries of the theory that have been broken by the regulators. 

For fixed total momentum _P+, a large number of particles must carry a very large light-front 
energy, since each particle is forced to carry a small positive Although this would also be true 
of massive particles in an equal-time quantisation, the effect is much more pronounced in the light- 
front formulation, because of the peculiar inverse relation between momentum and energy (|^). In 
particular, one notes that in the effective potential V, higher powers of the field Afr(x) would couple 
precisely to those components of the wavefunction with a large number of M-quanta. This is our 
physical justification for truncation of V to the shortest Wilson loops. Higher loops may be added 
as greater accuracy is required in the low-energy theory. 

There are potentially massless degrees of freedom associated with the continuum fields A± = 
(A" ± A^^^)/ \/2. The behaviour of these modes is subject to the gauge fixing and regulation of the 
fc+ = region, and is discussed in the next section. We will impose periodic boundary conditions in 
x~ [^Sj, with period C, taking the longitudinal momentum space continuum limit £ — s- cx) for fixed 
a. This gauge-invariant discretisation is also a useful starting point for the numerical investigation 
of the hght-front Hamiffonian ||l9| (DLCQ). 

3.2 Gauge Symmetry 

The theory @ possesses the fields AIr{x~^ , x~ A-i^{x'^ ,x~ and A-{x~^ , x~ ,x), together with 
the unitary gauge symmetry at each transverse lattice site x: 

A„(x) ^ C/(x)A„(x)C/t(x) + i(a„C/(x))C/t(x) (8) 
Mrix) C/t(x)M^(x)J7(x4-af) . (9) 
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The gauge fixing of this theory may be performed following similar treatments in Refs. pO| . With 
periodic boundary conditions on , we can pick the gauge = 0. This docs not generate a 

Fadeev-Popov Jacobian since it is (almost) an axial gauge. This leaves the zero modes 

A°_{x+,x)= [ dx-A^{x+,x-,x) , (10) 
which may be further gauge fixed to diagonal form 

A°_,,, iG{l,2,...,N} , J2a°_,, = 0, (11) 

i 

by -independent gauge transformations. It is convenient to define 

= ^A'-u ■ (12) 

The gauge fixing to diagonal form ( pT[ ) generates a Jacobian, which is just the invariant group 
measure for SU{N), proportional to 

[](sin(c,-c,)7r)2 (13) 

hi 

at each site x on the transverse lattice. This produces a repulsion between the eigenvalues A^ 
will turn out to be a constrained field which may be eliminated by its equation of motion in terms 
of the Mr{x.)'s and q's. 

There are residual global gauge transformations 



^m, -0, (14) 



where G {±2, ±4, ....}; Ci transforms as 

Ci ^ Ci — mi/2 (Gribov). (15) 

This generates Gribov copies. For pure gauge theory there are also the center global gauge trans- 
formations for each space direction x~ , x. For x" they are of the form 

U,, ^ 5,,e'2--""/^^ , i,j<N, (16) 

for n G {l,2,---,iV — 1}. For N — oo we can ignore the last matrix element and think of 
exp(27rx^s/£) • 1 for < s < 1. (More generally we will ignore the extra U{1) difference between 
SU{oo) and U{oo) when it does not matter.) Ci transforms as 

Ci ^ Ci — s (Centre). (18) 
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In the case that x^, say, is compactified with period Li, the center transformations are accordingly 
exp(27ria;^s/Li) • 1 for < s < 1. Under the latter, Mr(x) for f = (1,0) transforms as 



M^(x) ^ e2'^'''"/^iM^(x) . 
Thus a Polyakov loop (winding mode) W^(p.o)i of winding number p around Li, transforms as 

independent of ii. 

3.3 Confinement and Zero Modes 

The final gauge-fixed action is 



(19) 



(20) 



S — 5*1 + 5*2 + 5*3 



where 



Si 

S3 

J+(x) 



f dx+dx-^Re[Ti{d+Mr{x)D^M^^{x)}] -V^[M] 
j dx+dx- ^ Tr{(i?_^+(x))2} + Tr {A+(x)J+(x)} 
fdx+J2^{d^c,i^))' 

M^ix - af)D-Mr{x - af) - (D_Mr{x - af)^ Mr{x - af) 



Here D^cf) — d-(j) + i[A^, cf)]. The constraint equation of motion obtained from 52 is then 

2 



G2 



DiA+ = 



Formally inverting (Eq) one derives the longitudinal 'Coulomb' interaction at site x. 



S2 = J da;+rfa;-^^Tr{j+(x)(i?_)-2j+(x)} 



(21) 

(22) 
(23) 
(24) 

(25) 
(26) 
(27) 



We expect this to produce linear confinement in the longitudinal direction, because of the classical 
nature of this effect for the two-dimensional gauge theory at each lattice site. 
The degrees of freedom Ci are essentially winding modes around C, since 



Tr<^ P exp i / dx^ A^ 



(28) 
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in the present gauge. Therefore the classical linear confinement in two dimensions should push these 
modes to infinite energy in the £ — > oo limit. We can see this more formally in the large- iV limit. 
The Schrodinger representation of the kinetic term S3 (p4[), 

(29) 

evidently demonstrates that the modes q become classical at TV = 00, since G'^N is held finite. 
Introducing 

00 

4(x-,x)= i+(",x)e-2-'"-"/^, (30) 

n— — oo 



the longitudinal interaction (|27|) becomes^ 



G^£ y ,J . (31) 



^ (n-Q-|-Co)2 

n— — C30 ^ ^ 



The residual global gauge invariances, (|l5| ) and ([l^), mean that we can look for the classical minimum 
configuration in the fundamental domain < q < 1. From (|l|), we see that if J^jiO, x) 7^ 0, this term 
contributes an energy of order C. On the other hand, if J^{0,x) vanishes, say J^{n,K) ^ [n/Cy^ 
a finite energy is only obtained in the C 00 limit if s > 1/2. In this case, only n ^ 0{C) 
contributes to finite energy configurations and the Ci are negligibly small by comparison in the 
covariant derivative (States on which J^{0,x) vanishes also produce the necessary factor of N 



to make matrix elements of (|3l| ) finite in the large N limit; see Section 3^.) This argument treats 
the currents as essentially classical sources, and could be upset if quantum fiuctuations of A/j.(x) 
led to complete screening of the linear potential (|2^). Provided no condensation of the field AIr{x) 
takes place, as is assumed in the colour-dieletric regime, this cannot occur. 

The above arguments apparently justify use of the naive light-front gauge A_ = in the original 
work p. One possible caveat concerns the 0{£) contributions of q to the light- front energy, from 
the longitudinal Coulomb interactions (^ij) between link variables. Although these are infinite as 
C ~f 00, as demonstrated above, a principle value prescription operates between Coulomb self-energy 
and Coulomb exchange at zero momentum transfer, such that linear divergences cancel. It is 
plausible that the c,;, being relevant near the zero momentum region, might affect this cancellation 
in such a way as to alter the residual finite part. Very recently, some numerical evidence for this 



has been given for scalar SU{2) matter |21 (in the case of SU{2) there is only one value of i). 
In Appendix ^ we give an argument which shows that the only effect of the zero modes c; is to 
finitely renormalise the link-field mass /i^ via the singular Coulomb interaction (|3|). This is of no 



^The apparent divergence at n = 0, c; = Cj, can be avoided by using residual gauge freedom to gauge away 
^4-ii(^^) foi" each i at a specific light-front time x+, taken as the initial value surface = on which the constraint 
equations are solved |Q. 



12 



consequence for our approach, since /i is a renormalised parameter to be determined by physical 
measurements anyway. We will therefore proceed by dropping the variables Ci from now on. 

The charges J^(0,x) generate a;~-independent gauge transformations before gauge fixing, and 
the requirement that they vanish for finite energy states is the statement of colour confinement for 
the transverse lattice formulation at fixed a |^ . The mechanism for the confining linear potential in 
longitudinal and transverse directions is somewhat different in each case however. For the longitu- 
dinal direction it arises as the Coulomb potential for the 1 + 1-dimensional gauge field theory 
at each site x. In the transverse directions it arises through the Wilson-Kogut-Susskind argument 
for lattice gauge theories [|ll|, |2^. For colour singlet at each x, we must build a colour string with 
the (massive) A/r(x) link fields, the energy being proportional to the number of links in the string. 
Note however that the colour singlet requirement came from the continuum longitudinal dynamics. 
In later sections, we will show how one can verify these intuitive arguments by explicit calculation 
of the transverse and longitudinal string tensions. 

3.4 Energy and Momentum 

With the considerations of the last subsection in place, only Mr(x) remains as a physical degree of 
freedom. We derive the total light-front momentum and energy — J dx^ 6*+^, 

P+ = 2 /dx" ^ Tr{9_Af^(x)a_Afr(x)t} (32) 

P = |d:,-^y^[Af]-^(^Tr|j+(x)^J+(x)|-lTr{j+(x)}^Tr{j+(x)}) (33) 

with J^(x) now given by ( p5| ) with D- d-. The correction to the Coulomb term in ( |3^ ) is 
necessary at large- for SU{N) to ensure that the longitudinal gluon having no colour singlet 
part, does not couple to gauge singlet states. It will turn out to be important only for two-link 
states at non-zero transverse momentum P. 

Given the absence of modes Ci, we recover the residual local gauge symmetry generated by 
J+(0,x) 

Mr(x) -> C/t(x)Mr(x)[/(x + af) (34) 

with [/(x) independent of . Following the discussion of the previous subsection, the light-front 
energy ( |3^ ) is not finite unless the associated charge vanishes at each site x: J+(0,x) = 0. Thus, 
only states invariant under the residual gauge symmetry (^) will have finite energy. 
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3.5 Fock Space 



We impose canonical commutation relations at — 



Mr,ij{x ,X),{d^MsM{y ,y))^ = ^SilSjkS{x -y j^x.yf^r 



We work in longitudinal momentum space and transverse position space 



Mr{x+ = 0, a;",x) 



1 



dk 



aA,ij(fc,x),ap,fc;(fc,y) 



\/47r Jo Vk ^ 
= Sik Sji Sxp (5x,y S{k - k) 

= . 



(35) 

(36) 

(37) 
(38) 



Here, (oA.y )^ — {o.\)ji and A,/9 G {±1; . . . ± (_D — 2)} For clarity, we omit the + superscript on 
light-front momentum fc+ hereafter. The Fock vacuum |0) satisfies :P^: |0) = |0) = and is the 
physical vacuum. The Fock space operator a\{k,x) creates a mode with longitudinal momentum 
k on the link emanating from site x in transverse direction A. This Fock space is already diagonal 
in P+, and serves as a basis for finding the energy levels of . Since the glueball boundstates do 
not interact with one another at, N = oo, this in turn will diagonalise the invariant mass operator 
^ 2P+P- - P2 in the frame P = 0. 

The combinations singlet under (^) are the 'string' states which form closed flux loops on the 
transverse lattice, exempli gratia 



Tr|aLi(^i, X + al) a\{k2,x) a^j^(fc3, x + a2) 

■ a^_2{k4,x + a2 — al) a\{kQ,K. — al) a\{kQ,x)^ |0) 



(39) 



as illustrated in Fig. ^. The longitudinal momenta km are unconstrained except for km > and 
conservation, X]m=i ^™ ^ above example. 

4 Structure Functions. 

A standard gauge invariant definition in the continuum of the gluon distribution function G{x) in a 
hadron |^'(F+)) is || 

1 



G{x) 



dy-e-'-P^y {^{P+)\ F+^iO, y-,0) Oab F+^'{0, 0, 0) \^iP+)) (40) 



where F^" {x~^ , x^ , x) is the field strength tensor with colour index a,b,c£ {1, . . . , — 1} and 



C' = -pexpM^ dz-A+{0,z-,O)tc 



(41) 
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Figure 2: An example of a length 6 loop on the transverse lattice, showing also the longitudinal 
momentum carried by each link, corresponding to the Fock state (^9|). 
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Here, the tc are the SU{N) generators in the adjoint representation. In the hght-front gauge A_ = 
this definition (|4^) reduces to the expectation value of the gluon number operator. Note that the 
answer is in fact independent of the choice of x. Since the colour-dielectric formulation deals with 
a set of collective variables rather than with gluons, the transformation from one to the other being 
unknown explicitly, we cannot directly measure the gluon distribution. We recall however from (Q) 
and (H) that Z)_Afr(x) reduces to the field strength F^r in the continuum limit. It seems sensible 
therefore to make the substitution F^r D^Mri^) in the definition (pO[). We call the resulting 
gauge invariant distribution Gd to avoid confusion. In the continuum limit Gd reduces to G{x) but in 
general they are different. Since A/r(x) are lattice variables, Gd will evolve with the renormalisation 
scale defined by the edge of the Brillouin zone = vr/a. In the light-front gauge A_ = we have 

Gdix,A^) = —L_ J dy-e-'^''^y'{^iP+)\TT{d^Mrd^Ml}\^{P+)) (42) 
= (vl/(P+)|a^.(xP+,x)aA,.,(a;P+,x)|vl/(P+)) (43) 

Thus Gd has a simple interpretation: the expectation value of the number operator for link fields of 

longitudinal momentum fraction x. It trivially satisfies the momentum sum rule 

1 

dx xGd{x) = 1. (44) 







and the modes of Afr(x) may be interpreted as 'fat' partons with finite transverse extent. Bjorken 
scaling violations due to this finite extent will manifest themselves through the evolution with 



along the scaling trajectory. From the discussion of Section 2.1, one expects the scaling trajectory 
to pass from large to small values of parton mass /i^ as the lattice spacing a is reduced. Since a is a 
measure of the transverse size of the parton, the parton becomes fatter as it becomes heavier. As Aj^ 
is increased (a decreased) the glue distribution Gd is expected to evolve to smaller x. This follows 
simply from the fact that the small x region becomes enhanced the lighter the parton masses, since 
the free light-front kinetic energy (^ is ^ A^^/a; for a parton of momentum fraction x. We explicitly 
calculate structure functions later. 

5 The Large- Limit 

In pure gauge theory, the corrections to the large-A^ limit usually form an asymptotic series in 
1/N'^. For our problem, a dramatic simplification of the relevant Fock space takes place when one 
neglects these corrections. Firstly, in Ref. |^ and I it was shown that in the colour-dielectric regime, 
the Hilbert space of P = Wilson loops on the transverse lattice, together with the light-front 
Hamiltonian which acts upon them, undergo dimensional reduction; they may be replaced by an 
equivalent problem where all links in a given transverse direction are identified. One effectively 
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works on a transverse lattice with one link in each direction and periodic boundary conditions. This 
is a manifestation of Eguchi-Kawai reduction [ p5| in the light-front Hamiltonian framework. One 
may drop the argument x in all previous expressions, 



Mr{x ,x) ^ Mr{x ) , a\{k,x) ~* a\{k) 



(45) 



The proof of the dimensional reduction ( ^q ) does not encounter any phenomena analogous to center 
C/(l) breaking ||2^; however one must recall that the colour-dielectric regime, on which the proof is 
based, is only expected at large transverse lattice spacing. In this way, in D spacetime dimensions 
one ends up having to solve two-dimensional large- gauge theories coupled to 2^^^ complex adjoint 
scalars. The light-front quantisation of this kind of two-dimensional gauge theory has recently been 
developed in Refs. j22j, whose mathematical analysis we follow closely.^ 

The Hilbert space is further simplified by the fact that the light-front Hamiltonian propagates 
Wilson loops without splitting or joining at large- iV the amplitude for this is order Thus 
the light-front Hamiltonian may be studied in the space of connected Wilson loops, that is, states 



containing one colour trace, exempli gratia (39) 



The explicit light-front Hamiltonian (^3|) for the frame P = now becomes equivalently {D = 4) 
dx- ~^TtIj+^,P 

4: [ 0± 

- ^ Tt^mImIM2Mi + hermitian conjugatej + ^ /i^ TrjAfrM^} 

r 



E 



D-2 



N 



A 



3 



4A5 



Ti{MrM}MrMl} 
'Tr{MrM}}f + 



D-2 



N 



TT{MrMrM^M^} 



,D-2 



N 



E H 



Til Ml M.r Ml A'L, 



o-=±2,ct' = ±1 



TrjMiMtj TrjAfaAftj 



J+ = 



A ir d- Ml + Ml A ir 



(46) 
(47) 



In Section 7.1 we will show how P ^ may be treated in the same dimensionally reduced framework, 
by adding appropriate phase factors to matrix elements of ( ^ ) in the dimensionally reduced Fock 
basis. We have shown in ( p6| ) terms of the effective potential corresponding to Fig. |l]on the transverse 
lattice, with coefficients ,\-|^ _\2, A3, A4, A5, and /3. These parameters are to be fixed by a scaling 
analysis. Note that the multiple Trace operators coupling to A3 and A5 do contribute in leading 
order of N, but only on 2-parton Fock states. The dimensional reduction in the large- limit was 

■^The dimensional reduction discussed in those papers, which is an approximation, is not to be confused with the 
large-Af reduction discussed here, which is exact. 
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not discussed by the authors of Refs. j^, ^, but was used imphcitly in the calculations performed 
there. 

The fact, that the large-TV limit of the light-front Hamiltonian simplifies the calculations of 
physical observables in this way is to be contrasted with other quantisation schemes. A path integral 
evaluation of Green's functions 



S[A,] (48) 



should be dominated exactly at = oo by an x-independent classical saddle point configuration of 
the fields A^, the 'master field' |2^. This is motivated by the factorisation property 

( Tr{/l^j^p.-, • • • Af_i^} TrjA^^^j • • • ^/j„, }) 

= {Tv{A^,A^,---A^J) {Tv{A^^^,A,^^,---A^J) 

+ ^ { Tr{A,,,A^, ■ ■ ■ A^J Tr{A^„^,A,,„^, • • • A^J)^ (49) 



where (• • •)^ is the connected amplitude. Eqn. (49) can be understood by the usual planar diagram 



counting rules which attach a factor ]\[^-'^g-p diagram with G handles and P boundaries. 
The expectation value of two Wilson loops has a term of order N'^, from the disconnected disks {P — 1 
for each) covering each loop, plus a term of order N*^ from the connected amplitude, the cylinder 
connecting the two loops. Thus at = oo, where Eguchi-Kawai reduction in the lattice theory 
takes place, it appears that one can only calculate disconnected Green's functions and the vacuum 
energy. In an equal-time Hamiltonian approach, eigenvalues are also dominated by the vacuum 
contribution of order A^^ (amplitude A^ to create or annihilate a zero momentum Wilson loop from 
the vacuum). It is therefore not possible to extract the physical spectrum at A^ = cx), since this 
involves contributions down by order 1/A^^ compared to the vacuum contribution. Physical states 
interact with the vacuum by emitting or absorbing zero-momentum Wilson loops with amplitude of 
order A^'^ (amplitude 1/A^ to split or join and amplitude A^ to emit or absorb from the vacuum). 

It is then clear why the transverse lattice light-front Hamiltonian framework is uniquely suited 
for calculations of observables — as opposed to vacuum structure — in the large A^ limit. The triviality 
of the light-front vacuum implies that there is no order A^^ contribution. Non-trivial contributions 
begin at order A*^, corresponding to the free string- like propagation of a Wilson loop. In this way 
the spectrum of physical states is directly computable at A^ = cx3, where a dimensional reduction in 
the number of degrees of freedom takes place in the lattice theory. However, unlike the path integral 
vacuum problem where a reduction to zero dimensions takes place, in the light-front transverse 
lattice approach dimensional reduction takes place only in the transverse directions. The statement 
of large- A^ (Eguchi-Kawai) reduction in D spacetime dimensions for physical observables, rather than 
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vacuum structure, is that one can solve for observables in a field theory of D ~ 2 fewer dimensions 
in the large N limit. 



6 String Tensions 
6.1 Transverse Tension 

To measure the string tension aT in the x = {x^,x'^) directions, consider a transverse lattice with 
n = {Li/a, L2/a) hnks and periodic boundary conditions. We construct a basis of Polyakov loop 
winding modes in Fock space in the frame P = that wind once around this lattice in both the 
and x^ directions. The shortest such loop has length |n| in lattice units. One may extract from the 
lowest invariant mass eigenvalue the lattice string tension aaT via 

w a^crf, |n|2 (50) 



for sufficiently large |n|. The form (50) will be imposed as one of the scaling requirements, id est 
ax must not to vary with a or direction n/|n|. Because of large- iV reduction, this procedure is 
equivalent to using Polyakov loops of winding number n, on the transverse lattice with one link 
in each direction. The total winding number in each direction is conserved due to center symmetry 
(pO|). In the large- limit, we study single Polyakov loops for the same reason that splitting of 
Wilson loops W^(o.o) is suppressed by Thus n labels disconnected sectors of the Hilbert space. 

In sectors of non-zero winding |n| 7^ 0, there are further operators one could add to the effective 
potential T^[M] at the same level of approximation. They will in general be necessary for finiteness, 
as well as to achieve scaling in these sectors. In particular, when the Hamiltonian ( ^6|) is put into 
normal-order, two kinds of quadratic terms are generated, with divergent coefficients and c^, in 
the C —^ 00 limit: 

*Tr{4(fc)a,(fc)}; (51) 
r{4(fc)} TT{ax{k)} . (52) 

The first dH) is a link-field mass self-energy, and is dealt with by straightforward mass renor- 
malisation, leaving the finite part js^. The second ( ^^ is suppressed by 1/A^ on all states except 
Tr|a^(P^)|, the one-parton state of winding number A. To avoid the divergence of the mass of 
this state, one must add a new counterterm to given by 

-^Tr{MjTr{Aft} , (53) 

with an appropriate coefficient to cancel c^,. It is then possible to have a finite part left over, which 
will appear as a new parameter in the theory. In this case one loses one unit of predictive power 



'0 

Cw f . 
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from the measurement of energy levels in the winding number ±1 sector. Evidently, there are many 
similar multiple Trace operators one can add to the Hamiltonian, which act in this way on n 7^ 
sectors. They are quite analogous to the multiple trace operators, of which A3 is an example (^), 
which can appear in the n = sector (in that case winding number is conserved 'in each Trace'). 
However, ( |5^ ) is the only one required on the grounds of finiteness when the Hamiltonian ( p6| ) is put 
into normal order. In a winding n sector, multiple Trace operators that contribute in leading order 
of N must contain at least \n^\ powers of A/r(x). Since we try to extract ar from large |n|, we 
will neglect this possible proliferation of parameters to a first approximation, since we neglect higher 
powers of Mr(x) in the effective potential l/[Af]. 

6.2 Heavy Quark Potential 

In addition to using Polyakov loops of nonzero winding number in the transverse direction, we can 
measure the string tension by from the heavy quark potential [lO[ |. The advantage of this approach 
is that the heavy quark potential can be measured in the x^~^ as well as in the transverse directions. 
The main disadvantage is that this quantity is numerically difhcult to compute. 

For convenience, we will start with a heavy scalar field 4'{x^ , , x) of mass p in the fundamental 
representation of the colour group (one can just as well start with the Dirac equation and take the 
heavy limit, the result is the same). In addition to the pure glue 'link-link' interactions, the most 
general action to fourth order in the fields is 

J dx- dx+ {Do.(t))^ D"'(t) - p^(j)^(t) - ^^'"^3^ 0^ {M}Mr + MrMl) (54) 

where 

D^(t>'^ = dc,(t>^ ~ i<\>'^A^ . (55) 
The heavy field contributes to the gauge current J" in Eqn. ( |25| ) 

4" avy = (^5"'^) '^^ + *0 (^"'^)^ ; (56) 

note that (\)<^'^ is an iV x colour matrix. Since we are interested in the p 00 limit, we have not 
included any hopping terms like 0t(x)Afr(x)0(x + ar) and (f>{x.) is pinned to site x of the lattice. 
Also, we have not included any "pinching terms" like (f)^ (fxj)^ (f) which, for large iV, acts only on states 
consisting of two 0's. 

Now for some kinematics: Let P^^jj represent the full 2-monientuni of a system containing h 
heavy particles. It is convenient to split the full momentum into a "heavy" part plus a "residual" 
part P", 

Pf^ii = phv^ + P" , (57) 
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where is the covariant velocity of the heavy quarks, v"va = 1- The fuh invariant mass squared 
(P = 0) is 

= 2P+^^P^r^^ = {hpf + 2hpv+p- + 2P+ {P~ + hpy-) (58) 

Our choice of w+ is arbitrary and it is convenient to choose it such that P+ = 0. Consequently, 
u+P^ is just the shift of the full invariant mass A4 due to the interactions: 

M ^ hp + v+p- +Oil/p) . (59) 

Thus, u+P~ is the usual energy associated with the heavy quark potential; this will be our Hamil- 
tonian. 

We define creation-annihilation operators associated with the heavy field: 

dk 



(&(fc) e-*^"^""-*'^^ +dt(fc)e™-^°p+»fc^ ^ . (60) 



v47r \/ pv+ + k 

The creation/annihilation operators 6(fc) and d{k) have the usual commutation relations. The e*''""^" 
term removes an overall momentum pv" from the 2-momentum. The contribution of the heavy 
particle interactions to the "residual" momenta P" in Eqns. ( |3^ ) and (|3^ ) are, to leading order in 
large p and N, 

/oc 
dkk{bHk)b{k) + :d{k)d\k):) (61) 
-OG 

^h-eavy = -/ dk {bHk)b{k) + ■.d{k)dHky.) 

f dk,dk,dk,dkJ-^^^±^^—^^^^^bHk,)d\k2)dikMk3) 

An J (fci - ksy 

1 f dkidk2dk^dkA , , , , ^^ 

+S{ki - k2 -ks- ki) (Hi^s + Hs^i)] (62) 



with 



"^2^2 — 



G'(fc2 + fc3) r 
{k2-k^Y a^-^N ^ 

(■.d{k^)ax{k2)a\{k^)d\k^): + b\k^)a\{k2)ax{k^)b{ki)) (63) 

G'(fc2-fc3) T \ 

[k2 + k^Y aD-^Nj 

d{k^)a{{k2)a\{k^)d\k^): + b\k^)a{{k^)al^{k2)b{ki)) (64) 



^3-1 = ^Ls- (65) 

As an example, let us construct a state with two heavy sources with rest frame separation L in 
x^~^ (this corresponds to an x~ separation of 2v^L): 

|2) = [ dx- (j)\x+,x- +Lv-) (l){x+, x- - Lv-)\0)c-^'P^-^" (66) 
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1 r°° 

= -^J dkh'<{k)S{-k)e''^'' ^ . (67) 
The associated energy shift is 

(2|^+P-|2) G^iV|L| 
(2|2) = ■ 

Thus, as expected, we have hnear confinement with string tension of 1/4 in units of G'^N. 

Now, let us investigate the bound state equation for the general case. A state containing two 
heavy sources with rest frame separation L has the form (suppressing orientation indices A, p): 

IV) = ^-^==J^ dki---dkpi;{ku...,kp) J d/e^^^'" 

■b^l - K/2) a^ki) ■ ■ ■ a^kp) d^-l - K/2) (69) 



where P{l^y. = is the total longitudinal momenta of the link fields. The associated bound state 

equation has the form: 

p-v+i;{k,,...,kp) = :^V(fci,...) 



Stt J ^/k^ {k-i - k'Y 
G^Nv+ r dk[dk'^ k'^^k[ ,,,,,,, , -iLv-(k'+k') 



("link-link" interactions from Eqn. (46)) 

-h "source- link" interactions on the other end of the (70) 
string (complex conjugate of above). 

Note that the first term scales with momentum as {P-^^-^]^^ , the r terms scale as {P^^]^'^ , and all 
other terms scale as {P^^]^ ^ ■ Thus, the longitudinal momentum scale is determined as a balance 
between the first term and the remaining terms of the Hamiltonian. As opposed to the closed string, 
an instantaneous-annihilation process via can act in the p ~ 2 link sector and the "pinching" 
terms {exempli gratia A3) do not come into play; otherwise, there are no changes in the "link-link" 
interactions from the closed string case. 

There is a variety of possible cutoffs for ki that one might choose. The cutoff used here Pj^j^ < 
w+vG^iVPniax is convenient because the "link- link" interactions, which act on the interior of the 
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string, conserve P^l^-^]^- Holding the other momenta constant and ignoring particle number changing 
interactions, we see that fci has the limiting behaviors^ 

limV(fci,...) oc fcf , 0</3<l/2, -ii— =/3tan(7r/3) (71) 

* « i^-^j 

The first relation is exactly the same "endpoint behavior" that occurs for particles in the interior 
of the string. The second relation ( [72|) implies that the eigenvalues of v^P^ converge as 1/ (Pmax)^ 
for T — and 1/ (Pmax)'^ for r nonzero. 

We solve the bound state equation (It^) by rescaling and discretising the longitudinal momenta 



h ^ ~k, , h e {1/2, 3/2, . . .} , (73) 

where iiTmax > ^Unk ~ X)i ^4 maximum total discretised momentum. The ratio -Pmax/^max 

parametrises the coarseness of the discretization. The discretised bound state equation, expressed 



in terms of the dimensionless parameters Pmax and Ly/G^N, is linear in the couplings 



Vcm, -^=. -^=^ , (74) 

which all have units of mass. For large L, the wavefunction becomes oscillatory. In practice, this 
forces us to extract the heavy quark potential from relatively small longitudinal separations: of order 
a few transverse lattice spacings. 

7 Lorentz Invariance 

One of the advantages of light-front quantisation is that boosts — both longitudinal and transverse 
— become simple. However, our transverse lattice regulator, along with our use of effective degrees 
of freedom, break transverse Lorentz covariance. Instead, we will demand Lorentz covariance under 
transverse boosts for the lightest glueball states; we will use this to help determine the effective 
potential. 



In Section 7.1, we will construct states with nonzero transverse momentum and discuss the 



associated calculations of eigenstates of P .In Section 7.2, we will discuss how these boosted states 
can be used to help determine the effective potential. 

7.1 Non-zero transverse momentum P 

For the sake of clarity, we will show the transverse coordinates x for the link fields a\(k,x.). The 
transverse coordinates will later be dropped: we shall see that non-zero transverse momentum is 

^ The 13 here is conventional notation for the longitudinal endpoint behavior, not to be confused with the plaquette 
coupling in (^). 
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equivalent to introducing various phase factors in the Hamiltonian in the Eguchi-Kawai reduced 
theory. In the foUowing, we wiU consider vanishing winding number n = 0; the extension to nonzero 
winding is discussed in Appendix |^. 

A generic p-hnk state of transverse momentum P has the form: 

P)) - ^ E Tr{al^ (fci, XI + y) al (fc^, x^ + y) • • • 4^ (fc„ x, + y)} |0) . (75) 

where the winding number vanishes Y^^=i -^i — and y is summed over hnks of the transverse lattice. 
The transverse coordinates must obey the formula 

Xi = Xi_i + aAi_i , i<l<p. (76) 

Next, we must choose a phase convention for the overall phase x as a function of xi,X2, . . .. We 
have found a longitudinal momentum weighted center of mass 

to be especially convenient. For this choice of x, the component Af_r of the angular momentum- 
boost tensor which generates Lorentz boosts in direction r and obeys 

[p-,M^r] = iPr , = , [P\M-s\ = -iglP^ , (78) 

takes on a particularly simple form: 

e-'^"*^-a^(fc^x)e'^"*^- = a\{k,^) e''^^-'' (79) 
e"'^"*^-"- |*(P+,P)) = |*(P+,P -bP+)) . (80) 

In the computation of large- iV matrix elements of an operator O 

(^'(P+,P)|0|*(P+,P)) = (*'(P+,0)|e'^"*^-'-/-f'^C'e-'^''^^-'-/^^ l*(^'+,0)) , (81) 

one obtains phase factors in the Eguchi-Kawai reduced theory. These phase factors correspond to 
any shift in the weighted center of mass ([77|). For example, the interaction 

C = E Tr|a^^(fci,x) a^(fc2,x - a\) aj^(fc3,x - ap) a^p(fc4, x)| (5(fci + fc2 - ^3 - ^4) (82) 

X 

picks up the phase factor 

giP-M_./p+ g-iP-M_./p+ ^ ^i^_±M P • (a - p) ^ . (83) 

We can see from this example that the phase factors are functions only of longitudinal momenta and 
orientation indices A and p. Thus, Eguchi-Kawai reduction of the theory remains. 
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Although we have constructed states of definite transverse momentum, the construction of the 
operator P itself is somewhat problematic. Consider the most natural definition, 

/ dk a\{k , X + ar) — a\{k , X — af) a\{k,x.) , (84) 

which has the correct classical continuum limit. When we apply this operator to a closed colour 
loop, for example (^), we obtain a state that is no longer a closed colour loop and whose energy 
therefore diverges. Evidently, the correct form of P is quite complicated. 

In numerical work, it is advantageous to work with a real-valued representation of the Hamil- 
tonian. We note that, in the x'' direction, the boost operator breaks the transverse reflection 
symmetry : — s- —x^ of the theory. The boost changes sign under transverse reflections. 

For a boost in the direction, we introduce a basis with states that are even and odd under 7"". 
We use states of the form: 

\ijA) = cos f ) (1 + V'') |\E'(P+, 0)) (no sum over r) (86) 

\ijB) = sinf — (l--p'')|*(P+,0)) (no sum over r) (87) 

and the associated phase factors ( ^3|) become sines and cosines. The states |5'(P"'", 0)) are constructed 
to have definite charge conjugation C |*(P+, 0)) = ± |^'(F+, 0)). 

7.2 Dispersion Relations. 

We now indicate how a scaling analysis can be carried out, by enforcing 5*0(3, 1) Lorentz invariance 
of physical observables. There are a number of relativistic criterion which one might apply. With 
only the spectrum of glueball mass ratios, as a function of the couplings in V^[M], one might try to 
obtain Lorentz multiplets. This is a rather narrow test however, since it pays no attention to scalar 
boundstates. A better measure of relativistic covariance is the glueball dispersion formula. Perfectly 
relativistic glueballs should satisfy 

2P+P- =M'^ + IPP . 



Transverse lattice glueballs will in general satisfy 

2P+P- = G^N {Ml + Ml a^IPp + Ml a*|P|'' + • • •) . (89) 

We have used G^N as the overall scale, so that the measurable coefficients M^ above depend only 
upon the dimensionless ratios 

to2__m!_ i_ ^ ^_ I ^_ P /go) 
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In order that (89) agree with ( p8| ) one can tune the couphngs so that = and a?G'^NM.\ = 1 
for each glueball in the fow-lying spectrum. Due to the fact that there are transverse lattice hopping 
terms in P~, even when V = the dispersion relation automatically becomes asymptotically 
quadratic for small aP . For this reason, the condition M\ = is more difficult to use to fix V . The 
condition a^G^NAi\ — 1, which is more useful, ensures that the speed of light is isotropic. The 
rotational invariancc of the speed of light in the transverse directions was discussed in Ref. Q . In 
order to ensure that the speed of light in transverse directions ct equals that in the x^^^ direction 
cl (set to one by convention), one must determine the dimensionless combination a?G^N. 

One way is to measure both the longitudinal and transverse string tensions. This also allows 
one to express the glueball spectrum in terms of the string tension which is essential for making an 
absolute prediction of the glueball masses. Suppose we make the following three measurements of 
dimensionful quantities 

viG'^N ^ ca (91) 
V2G^N = a^cr^ (92) 
v^G^N ^ a (93) 

where wi, f2, are measured dimensionless numbers, c is an unknown constant, and a is the 
continuum string tension. Eqn. (|9l| ) comes from a glueball mass measurement. Eqn. ( |92|) comes 
from the transverse winding spectrum (Section |6.l|) . Eqn. (|93| ) comes from the asymptotic heavy- 
source potential in the x^^^ direction (Section 6.2). We have also demanded that the longitudinal 



and transverse tensions are equal ot — ^ a. Scaling means that the right hand sides of the 
above equations are not multiplied by dimensionless functions of a^G^N . We now show that this is 
enough to solve the theory completely. 

We use the slope of the dispersion relations at P = for a few low-lying glueballs to estimate 
the scaling trajectory of V , by demanding isotropy of the speed of light. From ( ^8|) and (|8^) one 
needs to determine G^Na? in order to compare transverse and longitudinal scales when comparing 
the speeds of light. From ( p2|) and (^) we can deduce 

G'^Na^ = V2/vl (94) 

There are a number of self-consistency checks one can perform, such as examining deviations from 
quadratic dispersion (non-zero Al2); the Lorentz multiplet structure in the glueball spectrum, and 
rotational invariancc of the heavy source potential. Similarly, we can now measure the lattice spacing 
a and glueball masses c in units of a. 
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8 Calculations 



We now implement the procedures described in the first half of the paper for the case of 2 + 1- 
dimensional pure large- A'^ gauge theory. The main interest firstly is to compute glueball masses 
in physical units. These can be compared with results from established techniques to ascertain 
whether the transverse lattice theory produces correct results, at least to the level of approximation 
we use. However, it is worth bearing in mind that, apart from confirming known data, the light-front 
Hamiltonian approach is uniquely suited to the calculation of hadronic wavefunctions in a general 
Lorentz frame, and it is the latter which chiefly motivates us. We will exhibit the structure of these 
wavefunctions also. We have also begun analogous calculations in 3 + 1 dimensions, which are briefly 
discussed in the conclusions. 

8.1 2+1 Dimensions 

Yang-Mills theory in 2+1 dimensions is superrenormalisable. Rather than a logarithmic Coulomb 
potential, a linearly confining interaction appears to be dynamically generated Therefore it is 
rather similar to 3 + 1-dimensional QCD. It provides a useful test-bed for non-perturbative calcula- 
tions, given that two space dimensions can be more efficiently handled numerically than three. The 
discussion of the first half of the paper can be trivially modified for 2 + 1 dimensions (the reader is 
referred to I for further detail) . Our light-front Hamiltonian ( ^ ) reduces in 2 + 1 dimensions to 

P- = J dx- Trj + /i^ TrjAfA-f^} 

+^ Tr{Af A/UO/t} + ^ Tt{MMM^M^} + (Tr{A/Af^})^ (95) 
J+ = iM d~ + d~ M (96) 

where now aG^ ^ as a ^ 0, with the standard 2+1 dimensional Yang-Mills coupling in 
the continuum limit.^ Otherwise the analysis is unmodified (in particular, the scaling analysis of 



Sec. 7.2). The link field M{x^), of course, now points in one transverse direction only. Under 
large- A'^ reduction it depends only upon , and is equivalent to two real adjoint matrix fields. We 
are now ready to diagonalise P^ in the Fock space of Af(x), as a function of m, hjh, h- G^N will 
be used to set the overall scale for dimensionful physical quantities. It will be determined in terms 
of a physical scale by measurements of the string tension. The glueball spectrum is classified by 
li/l''^'', for angular momentum J^, x^ reflection symmetry Vi, and charge conjugation C. The latter 
two are exact while first can only be estimated approximately on the transverse lattice. We follow 



^In I, we did not clearly distinguish between the a-dependent transverse lattice coupling G and the continuum 
coupling. Moreover, the normalisation of the continuum gauge coupling in that work was nonstandard. Replace 
by 2g^ and by rrP' jl to convert expressions in this paper to those appearing in I. 
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Figure 3: Convergence of the lowest glueball mass as a function of \/K for a 4-link truncation 
and — Xi — 0. The gray points are from conventional DLCQ and the dashed line is a fit 
to 1.1834 + 1.6395/v^ - 5.0731/K + A.3513/ K^^'^. The black points are from a calculation using 
improved matrix elements and are fit to 1.17752 + 2.44169/7^-4.68183/^^. Thus, given data points 
in the range 10 < K < 20, the improved matrix elements are necessary for a reliable extrapolation. 

the computational methods established in Rcfs. for adjoint scalar matter in two-dimensions, 
together with various improvements described in the Appendices of I and in the next section. 



Our basic numerical method is DLCQ |l9| , where we discretise the longitudinal momentum _P+ = 
2ttK/C into odd half-integers in the bound state equation. K measures the coarseness of the dis- 
cretisation and ii' — !■ oo is the longitudinal continuum. This approach has several advantages for 
numerical work: the resulting Hamiltonian matrix is sparse and matrix elements of the interaction 
are easy to compute. The main disadvantage of DLCQ is that numerical convergence tends to be 
slow. In Appendix C of I, we introduced an approach which eliminated the most severe convergence 
problems, the l/_fC^^ and l/\/X errors, while minimising the remaining l/K error associated with 
the instantaneous interactions. We have made some subsequent refinements to the method which 
are described in Appendix Because of these improvements, we are able to extrapolate to infinite 
K with some confidence. For instance, in Figure ^ we plot the eigenvalue of the lowest glueball, 
calculated using improved matrix elements, as a function of l/K. It is clear that a fit to cq -f ci/iiT 
should work quite well for K between 10 and 20. 

In additional to the harmonic resolution K, we truncate the Fock space in the number of links 
p. We calculate the convergence in p-truncation numerically by looking at the convergence of the 
lowest eigenvalue for winding number |n| = 1 and ( A', p-truncation) = (19/2,3), (25/2,3), (39/2,3), 



8.2 Methods 
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0.02 0.04 



Figure 4: Convergence of the |n| = heavy source potential with l/Xmax for Pmax = 4, 5, 8 and a 
2-Hnk truncation. The couphngs are from the m — 0.180255 datum of Table |^ and the longitudinal 
separation is L = A/VG^N = 1.750/^. The fit function is from Eqn. @. 



(19/2,5), (25/2,5), (23/2,7). Then we fit the eigenvalue to the function 

c„.(..|)e=...|. 



(<J7) 



It is satisfying that the results for Cp, typically of order -1, agree with the analytic estimate from 
Appendix B of I. We assume that, for sufficiently large p-truncation, Cp is a universal constant; it 
can be equally well applied to the spectra, nonzero winding |n|, and the heavy source potential. 

Using this result for Cp, we calculate spectra for various (ii', p-truncation) — (18/2,6), (18/2,8), 
(20/2,6), (20/2,8), (24/2,6), (32/2,6) and extrapolate to the continuum using the fmiction 



Co + cie 



K 



(98) 



We use the same extrapolation procedure for nonzero winding |n| ^ 0. In this case we use 
(|n|,X,p-truncation) = (2,20/2,4) (2,20/2,6), (2,24/2,4), (2,28/2,4), (3,21/2,5), (3,21/2,7), 
(3,23/2,5), (3,27/2,5), (4,20/2,6), (4,20/2,8), (4,22/2,6), (4,26/2,6). After extrapolating in K 
and p-truncation using Eqn. (|9^), we fit the lowest eigenvalue to the form {aTaln])"^ / (G^ N) + ci + 
C2/|np as discussed in I. 

The most tricky calculation is the heavy source potential. Here we have to extrapolate in Xmax, 
-Pmax, and p-truncation. In addition, we must choose reasonable values for the longitudinal separation 
L. For example, in Fig. ^ we see that the eigenvalues are well fit to the function 



v+P' 



1.64873- 



3.08839 



0.64118 



(99) 



(/'max)' V^max. 

with an RMS deviation of 0.0088. Although we would expect a the leading finite Xmax error to be 
of order (Pmax/^max)^, wc find the quadratic term to be much more important. Thus, we shall use 
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Figure 5: Convergence of the |n| = heavy source potential with l/Xmax for various Pmax and 
p-truncations as given by (101). The couphngs arc from the m = 0.180255 datum of Table |^ and 
L = 4:/V(PN = 1.750/V^. The fit function is from Eqn. ( |l02| ). 



the fitting function 



V+P- 



Co 



Cl 



{Pn- 



C2 



Pr 



Kn 



cae 



CpP 



(100) 



In subsequent calculations, we fit to ( |100| ) using the values 



p-truncation, P„ 



(0,32/2,2,4), (0,32/2,4,4), 



(0,32/2,2,5), (0,34/2,2,4), (0,44/2,2,5), (0,60/2,2,4), (1,19/2,3,4), 
(1,19/2,5,4), (1,19/2,3,5), (1,33/2,3,4), (1,33/2,3,5), (1,49/2,3,4) 



(101) 



Note that for the heavy source calculation, the conserved quantity n, analogous to winding number, 
is the number of links minus number of anti- links, id est the transverse separation of the sources in 
lattice units. Example data points are shown in Fig. |^ along with the associated fit function 

2 

(102) 



2 95576 / P 

1.54997 + . + 0.838171 ) + 0.773708 e'^"" 



In order to measure the heavy source potential and test its rotational invariance, we calculate the 



potential for |n| = and LV G^N = 3, 4, and 6. Then we fit to the functional form |29 

v+P 



= Co lVcpn + 



Cl 



C2 



(103) 



Next we calculate the potential for |n| = 1 and L\/(PN = 3; we then demand rotational invariance 
for this datum. 



9 Results 
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Table 1: The scaling trajectory which minimises the test of Lorentz covariance. The couplings 
t, and the overall scale G'^N/a are all obtained from the test, to is equivalent to the lattice spacing 
degree of freedom. 

9.1 First Principles 

To search for a scaling trajectory that restores Lorentz covariance, we applied a test for certain 
observables. First, we investigated which observables would be useful in this respect. The natural 
candidates are the ct/cl ratios deduced from glueball dispersion relations, the Lorentz multiplet 
degeneracy in the spectrum, and rotational invariance of the heavy source potential. 

In order to determine which glueballs it makes sense to include in the ct / cl part of the test, for 
each low-lying glueball we determined the li trajectory such that ct/cl — 1 for that glueball. In the 
low-lying spectrum, only the result for the first excited state 0^^ was markedly incompatible with 
the rest. In fact, it is consistently imaginary. We show below why this artifact is to be expected of 
our truncation of the effective potential to fourth order. Therefore, excluding this state, we included 
the ct/cl ratio of the lowest seven glueballs: 0++,2++,2 ^,0 ,0;7^,2^ ,2 . These are the 
lowest seven glueballs for a broad range of k at fixed m, and also the seven lowest in the ELMC 
simulations 

Another possible test of Lorentz covariance is the degeneracy of the parity doublets for the 2*+ 
and the 2^^. However we find that there is no point in our parameter space where the 2++ and the 
2 are degenerate and we do not include this pair in our test. For the 2^ and 2 , we have 
the opposite problem: the two levels are almost degenerate for a vary wide range of couplings. As 
a consequence, it does not provide a very useful criterion for finding the correct scaling trajectory 
(although it was included). 

In addition, we demand rotational invariance of the heavy source potential based on measure- 



ments at |n| =0 and 1. As described in subsection 7.2, the longitudinal string tension measurement 
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Figure 6: The trajectory which minimises the test of Lorentz covariance; see Table |l|. 
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Figure 7: The shaded hnes represent the heavy source potential as a function of separation 
V'cr + |n|2a2) for |n| = 0, 1, 2, 3, 4. The dashed hne is from Eqn. ( |l03| ) and the sohd hne 
is the heavy-heavy Coulomb potential from Eqn. (|68|). The eigenvalues were obtained by extrapo- 
lating in Kmax, ^niax, and p-truncation; the couplings are from the m = 0.180255 datum of Table 0. 



determines the overall scale G'^N/a and, as with our other measurements, we assume that there 
is some error associated with this. In principle, one might infer the overall scale from demanding 
consistent dispersion relations of the glueballs, but we found this method to give unreliable results. 

The coupling trajectory we find from the resulting test is shown in Table ^ and Figure |^. 
As previously found in I, the magnitude of I3 always comes out much larger than that of li and I2, 
essentially infinite. Here we are able to confirm that this is required by Lorentz covariance. We also 
elaborate on this below. Although there are some fluctuations in the coupling trajectory as m is 
varied, it is nevertheless remarkably close to the trajectory found in I on the basis of a fit to the 
ELMC spectrum. 

In determining the trajectory, we inferred roundness of the heavy source potential from mea- 
surements at |n| = and 1. If we plot the heavy source potential for more values of |n|, we observe 
that rotational invariance of the potential is maintained; see Fig |^. A potential that is "oval" — 
that is, rotationally invariant to within an overall ctt/cl — seems to be generic for the transverse 
lattice pO|] . Consequently, the crucial test is whether one obtains consistent longitudinal versus 
transverse scales (Tt = ctl = <J- 

Examining the glueball masses on this trajectory more closely, we find results for the 0"'""'' ground- 
state mass shown on Figure ^. It is to be anticipated that the becomes poor at large and small 
m. At large to, corresponding to large a (see Fig. pTl), the truncated effective potential cannot 
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Figure 8: The 0++ glueball mass and ^ values along the coupling trajectory of Table |l|. The size 
of the data points gives an indication of the error from extrapolation to the longitudinal continuum 
limit, id est the dominant error is from scaling violation on the transverse lattice. 
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Table 2: The glueball masses at m — 0.276 on the scaling trajectory of Table |i| at the minimum 
= 10.84. 

cope with the transverse discretisation errors. At small m, the wavefunctions are most sensitive to 
the the small fc+ region, and hence the longitudinal discretisation errors. It is also the edge of the 
colour-dielectric regime. The 0++ ct/cl ratio is 1 ± 0.15 all along the coupling trajectory, and the 
mass scales to some degree. Based on this, with a generous error for scaling violation, we estimate 
MQ++/y/a = 4.05 ± 0.3, in good agreement with the ELMC valueQ 4.08(7)(siai). 

At the best x^, we list the masses and ct/cl ratios of the other low- lying glueballs in Table 
The level ordering of the seven glueballs is in agreement with the ELMC result, and the numerical 
values are generally not far off. We note however that, with one exception, we cannot get the 
higher glueballs, along with the string tension and 0^^, to be especially rotationally invariant for 
our approximation to the effective potential. It is perhaps not coincidental that both the masses 
and ct/cl ratios of the higher glueballs come out consistently too low. The one exception, 2 '", is 
in good agreement with the ELMC mass 6.89(21) (stat). 

In Figures ^ and |l^, we plot the dispersion relation for the 0++ and 0^^ for a point on the 
scaling trajectory. This is a consistency check, to show whether optimizing the slope at |P| = 
leads to covariant behaviour throughout the Brillouin zone. Behavior in the figures is generic: the 
dispersion of the 0++ is extremely straight throughout the scaling trajectory and the dispersion of 
the 0~~, along with the other excited states, tend to follow a lattice dispersion cx \/l — cos(a|P|) 
more closely. 

We can also deduce the transverse lattice spacing a in physical units y/a, as one moves along the 
trajectory of Table |l] (see Fig. |ll|). This confirms that the lattice spacing gradually decreases with 

^The only other accurate value we know of comes from an SU(3) equal-time Hamiltonian lattice calculation [ po| . 
However, apart from corrections, it cannot be compared directly with our result since the authors do not measure 

the string tension, but give results in terms of a gauge coupling. 
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Figure 9: Dispersion relation for the 0++ up to the edge of the Brillouin zone. The sohd hne is the 
relativisticaUy correct dispersion relation and the dotted line is a lattice dispersion relation. The 
couplings are from the to = 0.180255 datum of Table |. 
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a|P| 

Figure 10: Dispersion relation for the up to the edge of the Brillouin zone. The solid line 
is the relativisticaUy correct dispersion relation and the dotted line is a lattice dispersion relation 
(normalised to the same slope at |P| = 0). The typically has a poor slope at |P| =0 and follows 
a lattice dispersion more closely. The couplings are from the m = 0.180255 datum of Table [|. 
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Figure 11: Variation of the lattice spacing a vs. to along the scaling trajectory of Table |l|. Also 
shown is a fit to 0.96 + 1.64to. 
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0.007 


0.889 
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0.76 




0.184 
0.716 
0.089 
0.68 
0.061 
0.871 



0.676 
0.189 
0.018 
0.118 
0.14 
0.002 



0.002 
0.009 
0.001 
0.023 
0.001 
0.005 



0.022 
0.027 
0.002 
0.053 
0.007 
0.061 



0.108 
0.059 
0.002 
0.056 
0.03 
0.061 



Table 3: The transverse structure of each glueball for K — 8, p < 6, a,t m — 0.134, li = —0.153, 
I2 — —0.105, ^3 = 100. Each column shows the probability for a glueball to contain a loop of given 
transverse shape. Information on the relative phases of different shapes can similarly be obtained 
from the wavefunctions. 
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Figure 12: The structure functions Gd{x) for K = 8, p < 6, at m = 0.134, Zi = -0.153, h = -0.105, 
Is = 100. This exhibits the longitudinal structure of each glueball. 
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m, and is always physically quite large in the colour-dielectric regime > 0. Further information 
can be gleaned from the glueball wavefunctions, which of course come for free in a Hamiltonian cal- 
culation. In Table | and Fig. |l| we exhibit the transverse and longitudinal structure of the glueballs 
at a representative point in coupling space. We see a characteristic rise of the G'd(a;)-distribution 
at small x, due to higher Fock states. However, one can show analytically that such distributions 
must vanish at x = for two-dimensional gauge theories [ pT[ (without Yukawa interactions), and 
therefore our colour dielectric regime rn^ > in particular. 

The coupling l^, according to the test for Lorentz covariance, is much larger than other couplings, 
which might appear unnatural at first sight. However the corresponding operator is also distinguished 
by the fact that at large N it couples only to link-anti-link Fock states (it 'pinches' the flux tube 
swept out light-front time x'^). More precisely, if the link -anti-link wavefunction is f{k^, — fc^), 
the ^3 operator has expectation value 



A complete set of link-anti-link wavefunctions / resulting from the longitudinal Coulomb interaction 
can be labelled by the number of zeroes of / |Q. The energy increases with the number of zeroes, and 
one might naively expect the 0^^ glueball to therefore be predominantly composed of a constant 
two-link wavefunction. However, this wavefunction and only this wavefunction makes a significant 
contribution to ( |104| ). A large therefore removes it from all low- lying glueballs, as is evidenced 
from Table |^ and fig. |l2|. An intuitive explanation for this follows from the fact that a wavefunction 
constant in fc+ implies tight binding in the longitudinal direction. On the other hand, the link- 
anti-link state has width a in the transverse direction. From Fig. |ll| we know this to be physically 
quite large (of order the glueball width!). Therefore it is perhaps not surprising that the Lorentz 
covariance conditions reject such an asymmetrical Fock wavefunction in low-lying glueballs. If this 
conclusion carries over to 3 + 1 dimensions, it will result in a quite different picture from the one 
suggested in Ref. Q (which did not include a I3 coupling). 

The absence of the link-anti-link state means that the 0++ groundstate glueball is predom- 
inantly formed from a linear combination of the 4-link operators Oi = TrjMtMM^M} and 
O2 = Tr{MMM^Aft|. In fact, we find that the (0++,0++) system is to a good approximation 
described by the symmetric and antisymmetric combinations {aOi+ j302, aOi— j302)- The operator 
O2 couples to ^2 < in the Hamiltonian and is the most significant hopping term on the transverse 
lattice. As a result, we find the symmetric combination lower in energy and with correct relativistic 
dispersion > 0. It then follows unavoidably that the higher-energy anti-symmetric combination 
will have an incorrect dispersion < 0. Since we tentatively identified the anti-symmetric com- 




(104) 
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bination with the 0^^, this explains why we always have a problem with this state. In fact, it is 
almost purely lattice artifact at present. It is likely that the effective potential must be extended to 
6-link operators before this problem can be resolved. 

9.2 Improving Higher Glueballs 

We have found in the previous section that, determining the glueball masses in terms of a scale 
set by the rotationally invariant string tcnsion(s), we can obtain a reasonable groundstate but 
cannot obtain Lorcntz covariant higher glueballs in general. We attribute this incompatibility to our 
truncation of the effective potential. It is natural to ask if one could find a coupling trajectory which 
improves the Lorentz covariance of higher glueballs at the expense of rotational invariance of the 
string tension. In particular, the longitudinal string tension is by far the most difficult measurement 
we make, and a potential source of systematic errors. To this end, we searched for such a trajectory, 
using the ELMC value of A/o++/\/ct as a phenomenological input to set the scale, cr was identified 
with ffT in this case. Although no longer quite a first principles calculation, the results are sufficiently 
impressive, we feel, to display here. 

Using a criterion which includes only the ct/cl ratio of the 0^"*" (whose mass is now fixed) 
and the , we found the trajectory displayed in Fig. In fact, along this trajectory not only are 
these two glueballs rotationally invariant, but a number of others also have a good ct/cl ratio. The 
trajectory is close to the previous one, but rather more stable. The scaling behaviour of glueball 
masses, shown in Fig. is also more stable. At the point of minimum on this trajectory, we find 
the spectrum displayed in Fig. In particular, it is worth noting that those glueballs whose ct/cl 
ratio is close to one have predicted masses in agreement with the ELMC result. This encourages 
us that, by adding further couplings to the effective potential to improve the Lorentz covariance of 
higher glueballs (from first principles), should bring their masses into line also. 

10 Conclusions. 

We have constructed an effective transverse lattice light-front Hamiltonian of pure gauge theory, ex- 
tending the discussion of quantisation and the large- iV limit addressed in the earlier works [Q, |[ |^ , 
and outlined a first-principles procedure for solving the eigenvalue problem. Not only does this 
formulation of gauge theory possess the usual appealing features associated with light-front quanti- 
sation — trivial vacuum, boost invariance, parton dynamics — but it also efficiently deals with the 
common problems found in this approach — zero modes, renormalisation, non-compactness/gauge 
invariance. It does this at the price of introducing an effective potential which must be tuned by 
demanding Lorentz covariance of eigenstates of the light-front Hamiltonian at large transverse lattice 
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Figure 13: The trajectory which minimises the test of Lorentz invariance involving only glueball 
dispersion formulae. 
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Figure 14: Scaling of the higher gluebaU masses M in units of -^/ct, as m is varied along the scaling 
trajectory of Fig. This gives an estimate of systematic errors for each glueball from truncation 
of V. 
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Figure 15: The low- lying 2 + 1-dimensional large N glueball spectrum obtained on the transverse 
lattice (o). Also shown for comparison, the ELMC simulations of Teper (+), extrapolated to iV = oo 
by a fit to A + B/N'^, with statistical errors The value of Mq++ / y/a was fit to the ELMC result, 
to set the scale. The other masses result from requiring rotationally invariant dispersion, quantified 
by the speed of light ratio ct/cl deduced from each glueball (in brackets). The estimated error from 
extrapolation in K and p is shown but an estimate of systematic error from truncation of V is not 
shown (see Fig. Oh. 
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spacing. Further remarkable simplifications occur in the large N limit. Only a 1 + 1 dimensional 
theory of connected Wilson loops needs to be solved. Equivalently, one solves a 1 + 1 dimensional 
continuum gauge theory coupled to complex adjoint matter fields which form freely propagating 
colour strings. 

We successfully tested the ideas by performing comprehensive calculations in 2 + 1 dimensions at 
large N , using an improvement scheme to extrapolate to the longitudinal continuum limit. We found 
agreement with conventional Euclidean lattice simulations js] , in those observables for which Lorentz 
covariance was attained. Furthermore, we could go beyond the path integral approach by obtaining 
the frame-invariant glueball wavefunctions. Given that the quantisation, regulators, elementary 
fields, and gauge fixing are different from those of Ref. the agreement is quite impressive. It non- 
perturbatively validates both the lattice colour-dielectric formulation and the light-front Hamiltonian 
quantisation of non-Abelian gauge theory at the quantitative level. This is therefore an important 
step in the application of light-front quantisation to hadronic physics in general. 

The obvious extension of large- A'^ pure glue calculations to 3 + 1 dimensions has already begun, 
and initial results are encouraging. For couplings in the light-front Hamiltonian ( p6| ) similar to those 
found for the 2 + 1 problem, we find 0++,2++, and 1^ niultiplets as the lightest glueballs (only 
some of the components of these multiplets were found in the original work |^). A search of the 
full parameter space for a scaling trajectory will yield estimates of the masses and wavefunctions. 
Given that 1/A^^ corrections are expected to be smaller than other systematic errors, the masses 
can be compared with various experimental glueball candidates, providing a useful alternative to 
the conventional ELMC estimates. 

Acknowledgements: SD thanks Prof. F. Lenz (Erlangen) and Mrs. van de Sande for hospitality 
during part of the work, and M. Teper for many helpful interactions. The work of SD was supported 
by a Particle Physics and Astronomy Research Council (U.K.) Advanced Fellowship and a CERN 
fellowship. 
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(a) (b) 
Figure 16: (a) Coulomb exchange graph; (b) self-energy graphs (contraction of (a). 



A A Zero Modes 

In this appendix, we argue that the only effect of the A_ zero modes q in the longitudinal Coulomb 
interaction (^) is to finitely renormalise the link field mass fi^. The two processes contributing to 
Coulomb exchange between link- variables ^/^(x) in light-front Hamiltonian bound state equations 
are shown diagrammatically in Fig. ^ (see for example Refs. for further discussion). Fig. p^(a) 
represents the instantaneous exchange of the non-zero mode, while Fig. |l^(b) is the self-energy 
contribution arising from normal-ordering the former. Neglecting Fig. p^(b) acts to regulate the 



singular behaviour of Fig. 16(a) at zero momentum transfer, by an equivalent principle value prescrip- 
tion Q. Let us denote the relevant Fock wavefunction of the two link variables by /(•••, a;i, X2, •• •) 
(or f{xi,X2) for short), for neighboring links in the colour Trace of a string state (e.gj39|), carrying 
momenta and X2P^ in the boundstate |\E'(P'^)) of total momentum P+. As explained in 



section 3.5, only neighboring links in the colour Trace have Coulomb exchange interactions in the 
large- limit. We can neglect labels for the link variables other than longitudinal momentum. Our 
basic assumption in the following is that in order to estimate the effect of zero modes Ci, which are 
constrained variables at large N whose distribution 0(ci, . . . , cn) is in principle known, we insert a 
c- number C/C wherever a quantity like Ci/ C would appear in the exact expression at finite C. In 
particular, C"s appearing in any expression should not be identified numerically with one another, 
but each represent something 0(1) in the continuum limit £ ^ oo. 

The contribution of the Coulomb interaction to the matrix element of the light- front Schrodinger 
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operator {Fock\2P+P^\^) is of the form^ 

r, N ni+m + C 



+ C){m + C) {ni - m + C)^ 
f{ni,n2){ni + m + C) f{m, rii + 712 — m){2n2 + ni - m + C) 



(105) 



[ yjirii + C)(m + C) V("-2 + C)(ni + 712 - m + C) 

Here Xi = ni/K, X2 ~ n2/K, and ni,n2,m ^ rii are non-zero integers less than integer K — 
£P+/27r. (If wc could set C = and take the limit C 00, we would recover the naive kernel of 
Ref. Pi, which employed the naive light-front gauge yl_ — 0.) For a finite free kinetic energy in the 



colour-dielectric regime, rii and 7^2 must be both 0{C). We Taylor expand (105) in C/C to obtain 

2ni + m! 1 



/(xi,x2) £ 



,t^„^ ^ni(ni-l-m') - m'Y 

/(ni,n2) / ' r^f, 

' zni + m —CI — 



v/?ii(ni + to') V V ni{ni+m') 

/(ni+TO',n2-TO') ^2n2-m'-C('l ^^"^ " 



^712 (n2 - m') V V '^2('^2 - m') 

+ 0{^) . (106) 



The higher corrections 0{C/C) can be dropped in the contimmm limit. Wc consider first the term 



at 0(C) in the curly brackets of (106). As £ ^ 00 it can only affect the finite part of / through the 
region of infinitesimal momentum transfer m' = to — rii ^ 0{1). In this region it simplifies to 

Cf{n,,n2)(- + -) 77^ 7^ (10^) 

\ni 712/ (C — m'Y 

^ ^ m' ——ni 

which, when folded with the appropriate distribution 0(ci, • • • ,cjv)^, is equivalent to shifting the 
mass in boundstate equations. The key point is the absence of cross-terms like 1/^x1X2- To 
actually calculate the mass shift would require 0, which we have not tried to determine in detail. 
In the case of SU{2) scalar matter, it should be possible numerically to check whether the result of 
Ref. 1^ can be reinterpreted as a mass shift. Note that our argument does not lead to any mass 
shift for fcrmionic matter. 



Turning to the term 0(C*') in the curly brackets of (106), this reduces to the continuum kernel 
of Ref. 1^ but for the contributions in the region to' ~ 0(1). To discuss the latter, we must specify 
more carefully the C"s. The correct expression for this region is (c.f. (^l|)) 

/ dc,|0(ci,---,cjv)p4A/(ni,n2)£ V 7 -^ + 0(l/£) (108) 

Jo „,^_„^ ic, -Cj-m) 



*Here we neglect to show mass ronormalisations which are known aheady from the naive analysis 
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where we have now exphcitly integrated over the zero mode distribution (f>, and used 
/(ni, na) - /(tt-i + to', - to') = to' A/(ni, ria) + 0(l/£^) 



(109) 



with A/(ni,n2) ~ 0(l/£). At first sight it appears as if the zero modes Ci spoil the cancellation 
of (|lO^ ) between positive and negative m' which would occur in their absence. However, using the 
fact that (j) has definite symmetry under interchange of any two c^'s, ones sees that after integration 
over Ci's, the cancellation does hold. Another, more physical, way to see why this cancellation must 
occur, is that a derivative term df{x)/dx would result in the continmmr bound state equation if it 
did not. Such a term would violate x^^^ — > —x^^^ symmetry in general. 

B Nonzero winding and nonzero transverse momentum 



The analysis of nonzero transverse momentum in Section 7.1 can be extended to nonzero winding 
number n ^ 0. Thus one can also study the Lorentz boost properties of the winding modes (we 
have not yet attempted systematic study of this). 

Let us start with a p-link state that winds once around a lattice with n transverse links. A 
generic P = state has the form 



y 



•4^(fcp,xp + y)}|0) (110) 



where y is summed over aWV — ■ ■ ■ ^ sites of the lattice and J^i = ^- By periodicity of the 
lattice we identify a|(fc,x) with a|(fc,x + an). The transverse coordinates are given by the formula 

Xi = Xi_i + a Ai_i , l<i<p. (Ill) 

Note that the number of nonzero contractions w ~ (^'(P+, 0)|5'(P+, 0)) must divide evenly. 



rf /w e -E, for all r. Now, we will boost this state (IIC) to some finite transverse momentum P: 



|*(P+,P)) = e*^"*^-/^^ |*(P+,0)) 

= ^^E^*''''^Tr{"^^(^i'y)"^^('^^'^2+y)'--"l^^^'^P+y)}l'^^ -^^^2) 

For nonzero winding, the phase condition introduced earlier ( [ttI ) has no particular advantage; instead 
we have chosen the convention Xi = x = 0. 

Since we have a periodic lattice, P can only take on discrete values 

aP'n'' 

r = — e Z (no sum over r). (113) 

27r 
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We find the inner product 

1^ , otherwise ^ ' 

Thus, as we change the transverse momentum, the number of states in a basis can change. Unhke 
the n = case, the boost operator is no longer self-adjoint due to the finite transverse volume of 
the system. 

In practical calculations, we prefer to work with real numbers. Thus, we introduce orientation 
reversals 

O Tr|at ^ (fci, xi)a^^ (fca, xa) • • • at,^ (fcp, Xp)| |0) 

= Tr{al^(A:p,yi)4^_^(fcp_i,y2)---<(fci>yp)}|0> • (US) 
We perform no cyclic permutations on the colour trace. We find 

yj = n-Xp_,+i (116) 
Consequently, Oe*^ O = e~'^ where B = P^M-r/P^ ■ It is useful to introduce the two states: 

|A(P+,P)) cos(B)(l+C')|*(P+,0)) 

|B(P+,P)) = sin(B)(l-O)|*(P+,0)) . (117) 

For states without symmetry under orientation reversal (^(P+, 0)|O|^(P+, 0)) = 0, we find the 
inner products 

(A(P+,P)|A(P+,P)) = (P(P+,P)|P(P+,P)) = 1^^' ^^j^L^^ (118) 
(A(P+,P)|i?(P+,P)) = (B(P+,P)|^(P+,P)) = (119) 

For states which are symmetric under orientation reversal (^(P+, 0)|O|^(P+, 0)) ^ 0, |j4(P+,P)) 
and |P(P+,P)) are equivalent up to a phase factor. In this case, we include only the |A(P+,P)) 
state in the basis. Occasionally, (A(P+, P)|A(P+, P)) = in which case we include |P(P+,P)) in 
the basis instead. 

C Improved matrix elements for DLCQ 

In this section, we will describe some refinements to the basic technique introduced in Appendix C 
of I. Our method is the following: when calculating the matrix elements of the Hamiltonian in a 
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many-particle calculation, we will use ordinary DLCQ commutation relations to calculate contrac- 
tions associated with spectator particles. To calculate the interaction itself, we take the DLCQ 
momenta of the interacting particles, project onto a smooth wavefunction basis (continuous longi- 
tudinal momentum), and calculate matrix elements of the interaction in this smooth wavefunction 
basis. Our modification of DLCQ essentially eliminates the \/^/K and errors associated with 

small momenta. The remaining errors include a (small) 1/K error from the instantaneous interac- 
tion l/(fc+ — p+)^, errors from DLCQ itself, and a 1/K error from production of particles 
at small momentum. In the following, we discuss improvements to the (2 particles) (2 particles) 
interactions, the (1 particle) — > (3 particles) interactions, and the heavy source interactions. The 
'particles' in the present context arc the link fields of course. 

C.l (2 particles) (2 particles) interactions 

Consider a typical 4-point interaction 

ki ^^^^ fcs 



K' K 



k2 -""^ ^ fc4 

with DLCQ momenta ki S {1/2, 3/2, . . .} and K — ki + k2 — + k^. In particular, the mass term 
is treated as if it were a 4-particle interaction: 

V{x,y) = ^(- + -^)s{x-y). (121) 



1 

We start by calculating eigenfunctions of the two-particle bound state equation 

M'^f{x,l- x) 
1 



dyH{x,y)f{y,l-y) 

ti^f{x,l-x) G^Nf'^^ f{y,l-y) / {x + y){2~x~y) ^ ^^^2) 



2^(1 - 47r Jo y/x{l-x)y{l-y) \ (x - yf 

where Ai + A2/2 + A3 = —G^Na^^^. Then we write the incoming and outgoing DLCQ states 
as linear combinations of the first K eigenfunctions and evaluate matrix elements of the various 
interactions V{x^y) in the eigenfunction basis. 

We introduce a basis of even and odd wavefunctions (labeled by subscripts s and s'): 

Ux) = {x{l-x)f[ co«(-^(^-f)) ' , (123) 

^ ^ ' ^ ^ " \ sm(7rs(a;- i)) , s G {1, 3, . . . , 25 - 1} ' ^ ' 

with (3 is given by = tan(7r/3), < /3 < 1/2. We define, for both even and odd wave- 

funtions, matrix elements of the various interaction terms Vg.s' = ^ dx dy (j)s{x)V {x^y) (j)si{y), 
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the inner product matrix E^^s' — J dx dy (l)s{x) (l)s' {y), and the two-particle Haniiltonian = 
J dx dy (j)s{x) H (x,y) (j)s'{y). In addition, we define the matrix Ti^s which maps the wavefunction 
basis onto the DLCQ basis 



T — 



kik2 



cos TTS 



sm 



(- (^)) 



■ 2K . 
fcl-fc2^ 



(124) 



where the subscript i E {1, 2, . . . K} denotes a two particle DLCQ basis state with momenta ki — 
i - 1/2, k2 = K-ki. 

This choice of basis wavefunctions has several advantages over our previous choice. First, this 
basis is numerically 'well-conditioned' while the previous choice was numerically unstable. Second, 
the highly excited eigenfunctions of ( |122| ) are themselves nearly sines and cosines. We now generate 
the improved matrix elements at the beginning of a calculation, store them in memory, and use them 
as needed. 

One calculates the various integrals using the formula 

dx x^'il - xY cos (2a {x ~ 1/2)) = ^fliV{pL + 1) /^+i/2(2a) (125) 
where /^(x) is defined in terms of a Bessel function: 



Thus, for example, 



3,{zl2)lz^ 



V/mC^;) = ./p-i(a;) +a^VM-Hi(a;) 



lim /^(x) 



r(^ + i)22m 



(126) 
(127) 
(128) 
(129) 



The greatest difficulty comes from the instantaneous interaction 



Stt 



dx dy 



{x + y){2 - X - y) (f>six) (y) 



(130) 



(131) 



y/xy{l - x){l - y){x - yY 

We cannot solve this integral analytically, but we can subtract the pole part aX x = y and solve the 
remaining finite integral numerically. Thus, we add and subtract the quantity 



Stt 



dx dy 
'I 



{x + y){2 ~x-y) {x{l - x)y{\ - y)f-^'^ 



{x - yf 



■ cos 



-{s - s'){x + y - I) 



±-cos 



-Xs + s'){x + y - I) 



(132) 
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which, using the identity 



(x - y)' 



= — hm 



dz z cos{z(x — y)) e 



(133) 



equals 



[fp{v-) [(1 + 4/5) fp{v-) - 2 
±f0{v+) [(1 + 4/3) f0{v+) - 2 } 



(134) 



where = 7r(s±s')/2. The numerical integration of (131) minus (132) is performed using the GauB- 
Jacobi method | |3^ which correctly handles the endpoint behavior of the integrand and converges 
exponentially fast. 

We solve the generalised eigenvalue problem Hxt = htExt where the eigenvectors xt have nor- 
malisation Xt'^Xt = St,t'- Using the lowest K eigenvectors, we define an eigenstate basis (denoted 
by subscripts t and t'), 

T^,t = T,,s ixt)s ■ (135) 

To make the transformation between the DLCQ and eigenfunction bases orthonormal, we perform 
a standard QR-factorisation [ ^2[ , Ti^t — Qi.t'Rt'.t, and discard the triangular matrix R (diagonal 
elements of R must be positive). Finally, we define the improved matrix elements by 



C.2 (1 particle) (3 particles) interactions 



(136) 



From the general bound state equations of two-dimensional gauge theories with pair production, one 
can show that as the momenta of any two particles go to zero the wavefunction becomes constant 
. Thus we define a basis of wavefunctions 

(137) 



{xy(\-x-y))' 



(pk.MMi^^y) = Ck.MM \ {{x + y)(l - a;)(l - y)f^^ 

, 



otherwise 



with normalisation 



1= dx 



l-x 



dyi 



'ki,k2,k: 



{x,y) 



(138) 



for three outgoing particles with longitudinal momentum fractions x, y, and l—x—y in the continuous 
wavefunction basis and ki e {1/2, 3/2, . . .}, K ~ ki + k2 + k^, in the DLCQ basis. We define the 
polygonal region Rki,k2,k3 surrounding {ki/K,k2/K) as shown in Fig. |l^. 
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K=3/2 




0.2 0.4 0.6 0. 



K=5/2 




0.2 0.4 . 6 



K=ll/2 




0.2 0.4 0.6 0. 



Figure 17: Region RkiMM ^'^^ ^ ~ 3/2, K = 5/2, and K = 11/2. The dots represent the points 
{ki/K,k2/K) and RkiMM the surrounding polygonal region. 

Consider a typical 4-point interaction 



V 



(139) 



We calculate improved matrix elements associated with interaction V{x, y,l — x — y) using the 
expression 

Vfci,fc2,fc3 = / dx dy(f>kuk2MiX'yWiX'y^'^-x-y) ■ (i^o) 
c/o Jo 

The integrals are calculated numerically, paying careful attention to the small momentum regions. 

C.3 Heavy source interactions 

Improvement of the matrix elements for the interactions between the heavy sources and the link fields 
is particularly challenging because of the oscillatory nature of the wavefunction for large separation 
L and the presence of the cutoff Pmax- So far, we have concentrated on removing the leading errors 
which are, as before, associated with small longitudinal momenta and arc of order l/(Jirmax)^^ and 
1 / V^max • At this point it makes no sense to include the finite L phase factors in the improvement 
scheme because we will also neglect oscillation of the wavefunction itself in the improved matrix 
elements. The leading error associated with the improved matrix elements is of order l/iiTmax- 
For (source-link) (source-link) interactions, we define wavefunctions 

"'^ fc- 1/2 < a; < A; + 1/2 



(j)k{x) = Ck 



X' 





otherwise 



(141) 



where k,k' E {1/2,3/2,...}, restricted by i^max > -K'link, is the DLCQ momentum of the link 
field neighboring the heavy source and x is its continuum longitudinal momentum rescaled by 
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Kaiavi/ {v^ Pmax) ■ The normalisation is 1 = J dx(f>k{x)'^. We calculate improved matrix ele- 
ments for the various interactions using 



Vk,k' = J dxdy Vix, y) 4>k' (y) (142) 

with the finite L phase factor replaced by its DLCQ value 



exp 777^ ^ exp — . (143) 

In the case of the mass term for the link field adjacent to the heavy source, we assign half of 
the contribution to the link-link interactions acting on the interior of the string, consistent with 



Eqn. (121), and the other half to the heavy-link interaction. Thus, we use 



V{x,y) ^ —^^ ^^"^ 5{x-y). (144) 

For the instantaneous A-^- exchange interaction, we must be careful about how we handle the 
l/{x — y)'^ pole and the i^max > -f^iink cutoff. For the off-diagonal part, we use 

^'^^---Mfc^mw^^P \jdxdy-=M^)^Ay), k^k . 



(145) 

If the l/{x — j/)^ pole were left inside the integral, there would be a logarithmic divergence when 
k — k' ±1. The diagonal part is obtained by adding and subtracting the 'self- inertia' 

UJ "^^^io ''VTyik-yf W 

where K' — Xmax — K\in\i + k. The integral is solved numerically with an implicit principle value 
prescription for the pole at y = k. This method for handling the self-inertias is discussed in Ref. 

Turning our attention to the (source) (source + 2 links) interactions, we proceed in an anal- 
ogous manner. Since the wavefunction becomes constant as any two momenta go to zero, we define 
wavefunctions 

( {xyY fci - 1/2 < X < A:i + 1/2 and 

c^k^M{x.y)^Ck,M\ {x + yyP' fc2-l/2<y</c2 + l/2 (147) 
[ , otherwise 

where fci, /c2 G {1/2, 3/2, . . .} are the DLCQ momenta of the two outgoing link fields created adjacent 
to the source. The normalisation is 1 = ^ dxdy (t>ki.k2 i^i 2/)^- We calculate improved matrix elements 
for the various interactions using 



Vkuk2= J dxdy(j)ki^k2{x,y)V{x,y). (148) 
As before, the finite L phase factors in V{x, y) are replaced by their DLCQ values. 
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